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Abstract 

The ring of Fermat reals is an extension of the real field containing 
nilpotent infinitesimals, and represents an alternative to Synthetic Differ- 
ential Geometry in classical logic. In the present paper, our first aim is to 
study this ring from using standard topological and algebraic structures. 
We present the Fermat topology, generated by a complete pseudo-metric, 
and the omega topology, generated by a complete metric. The first one is 
closely related to the differentiation of (non standard) smooth functions 
defined on open sets of Fermat reals. The second one is connected to the 
differentiation of smooth functions defined on infinitesimal sets. Subse- 
quently, we prove that every (proper) ideal is a set of infinitesimals whose 
order is less than or equal to some real number. Finally, we define and 
study roots of infinitesimals. A computer implementation as well as an 
application to infinitesimal Taylor formulas with fractional derivatives are 
presented. 
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1 Introduction 

If mathematics is the language of nature, one can imagine that the more re- 
sults are discoverable and describable using a given part of mathematics, the 
more faithfully that language will correspond to some given part of nature. 
We can hence imagine that there is a sort of weak isomorphism between that 
language and the corresponding part of nature. Therefore, if two different lan- 
guages are able to describe faithfully the same part of nature, we can also 
think that these two languages are, in some way, isomorphic to each other. So, 
because we are able to use mathematical analysis and actual infinitesimals as 

'University of Vienna, Faculty of Mathematics, Nordbergstr. 15, A-1090, Austria. Sup- 
ported by a L. Meitner FWF (Austria) grant (M1247-N13). 

t University of Vienna, Faculty of Mathematics, Nordbergstr. 15, A-1090, Austria. Sup- 
ported by FWF research grants Y237 and P20525 



f 



languages to describe nature, we can imagine that using the first one it would 
be possible to obtain a rigorous and modern model of the informal use of in- 
finitesimals. If this idea is in some way correct, this should actually be feasible 
without any non trivial background of mathematical logic. The theory of Fer- 
mat reals represents a possible formalization of this philosophical idea. Other 
possible approaches following this line of thought are: Weil functors [T71 ITS] , 
Levi-Civita fields I2r?] . Surreal numbers [5], geometries over a general 

base ring [2J, or the ring of Colombeau generalized numbers [SHHH]. Classical 
approaches requiring a non trivial background of mathematical logic are Non- 
standard Analysis and Synthetic Differential Geometry [T5I HZT] . In case 
the above philosophical idea sounds natural to the reader, an open problem is 
to understand, from a mathematical, social or historical point of view, why the 
latter theories, i.e. those requiring a non trivial background of mathematical 
logic, seem more powerful than the former ones. 

The ring *R of Fermat reals can be defined and studied using elementary 
calculus only ([IS])- It extends the field R of real numbers and contains nilpo- 
tent infinitesimals, i.e. h G *R^o such that h n = for some n € N>i. The 
methodological thread followed in the development of the theory of Fermat re- 
als has always been guided by the necessity to obtain a good dialectic between 
formal properties and their informal interpretations. Indeed, the ring is totally 
ordered and geometrically representable ( [HI HH [S] ) , to cite some examples. 

Every Fermat real x £ *R can be written, in a unique way, as 

N 

x = °x + y^ j a i - dt ai> (1.1) 

i=l 

where °x, a,, ai € R are standard reals, ai > a,2 > ■ ■ ■ > > 1, 04 ^ 0, and 
where dt a verifies the following properties 



dt a ■ dtb = dt ab 
(dt a f = dta M P g E>i (1.2) 
dt a = VaeM<i. 

The expression (|l.ip is called the decomposition of x, and the real number °x 
its standard part. The number a% =: u(x) is called the order of x and represents 
the greatest infinitesimal appearing in its decomposition. In case x G R, i.e. 
x = °x, we set tu(x) = 0. We will also use the notations uii{x) := and 
°Xi := at for the i-th order and the i-th standard part of x; Ui{x) := if x G R. 
The order uj(— ) has the following natural properties 



uj(x + y) = max [uj(x),uj(y)] 
1 1 1 



uj(x-y) u)(x) uj{y)' 
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whenever x, y are infinitesimals such that x + y =/= or x ■ y 0, respectively. 

Directly from (jl.ip it is not hard to prove that if k G N>i, then x k = iff 
uj(x) < k. For a G R>o U {oo}, the ideal 

D a := {x G *K | °x = 0, u(x) < a + 1} 

plays a fundamental role in fc-th order Taylor formulas with nilpotent increments 
(so that the remainder is zero). Indeed, for k G N>i we have that Dk = 
{x G *R | x k+1 = 0}, and any ordinary smooth function / : A — > R defined on 
an open set A of R d can be extended to the set 



•A = {a; G 'R d \ °x G ^} 



obviously obtaining a true extension, i.e. the same values at x G A. The men- 
tioned Taylor formula is therefore 

VheDi: f(x + h) =53^.^( a)> 

Ul<fc 

where a; G A is a standard point, and Df = Dk x . . A . . xDk- 

It may seem difficult to work in a ring with zero divisors, but the following 
properties permit to deal effectively with products of nilpotent infinitesimals 
(typically appearing in multidimensional Taylor formulas) and with cancellation 
laws: 

n 

x is invertible -4=> °x =/= 

(If x ■ r = x ■ s in *R, where r, s G R and i^O) =>■ r = s. 

Finally, the ring *R is totally ordered, and the order relation can be effectively 
decided, once again, starting from the decompositions: let x, y G "R; if °x ^ °y, 
then 

x < y <=> °x < °y. (1.3) 

Otherwise, if °a; = °y, then 

1. If u){x) > uj{y), then x > y iff °x 1 > 0, i.e. iff x > (from ([P]0 

2. If w(x) = to(y), then 

"xi > °yi => x > y 
°xi < °yi =^ x < y. 
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For example, d<3 — 3dt > dt > 0, and dt a < di& if a < b. This quick 
summary of some algebraic and order properties of the ring *R of Fermat reals 
can be considered as a first step toward its axiomatic description. A very simple 
model can already be guessed from the properties and (|1.2p . Indeed, we 
first introduce the ring R D [i] of little-oh polynomials, i.e. functions x : R>o — > 
R that can be written as x(t) = r + ^2f =1 cti -t ai +o(t), as t — > + , where A 6 N, 
r, Ox, • • • , &N £K and ax, . . . , aw € R>o- Then, in the ring R [i], we define the 
equivalence relation x ~ y iff x(t) = y(t) + o(i), for t —> + , and *R := R D [i]/ ~ 
is the related quotient set. 

For a more complete presentation and study of the model, see P2] HH [H| . 

In [10], the Fermat-Reyes theorem, which is essential for the development 
of differential calculus on *R, is presented. The Fermat-Reyes theorem states 
the existence and uniqueness of the smooth incremental ratio of every smooth 
function f : *U — > *R, that is existence and uniqueness of a function r : *U — > 
•R satisfying 

f(x + h) = f(x) + h-r(x,h) in *R V(x,h) £*U, 

where *U := h) G *R 2 | [x, a; + h] C *t/| is called the thickening of *[/ ([TJ 

[TU]). Here, the function / is more general than the extension from U C R to 
*[/ C *R of an ordinary smooth function defined on U and with values in R. 
The function / is a non standard smooth (or, more simply, smooth) function, 
i.e., by definition, / can locally be written as 

f(x)='a(p,x) Vx€'Vn'U, 

where a € C°°(W x V, R*) is an ordinary smooth function defined on an open set 
of R w x R v and p G *W is a w-dimensional parameter. The mentioned topology 
is that generated by extension of open sets, i.e. by sets of the form *U . For 
example, f(x) = dt + x is (non standard) smooth, but it is not the extension 
of an ordinary smooth function because /(0) = dt R, whereas any extension 
takes R to itself. The Fermat-Reyes theorem is well framed in the cartesian 
closed category *C°° of Fermat spaces, as this permits to develop a notion of 
smooth space and smooth function including also infinite dimensional spaces, 
e.g. function spaces like Man(A/, N) or integral and differential operators. The 
definition of °C°° is essentially a generalization of the notion of diffeological 
space ([2]), whose category C°° is the domain of the Fermat functor *(— ) : 
C°° — > m C°°. This functor generalizes the construction R n- *R and displays 
very good preservation properties, closely related to intuitionistic logic. For 
more details see [TUl • 

The main aim of the present paper is to study the ring of Fermat reals us- 
ing standard topological and algebraic structures. We will analyze interesting 
metric structures deeply related to the development of smooth calculus on *R. 
We will characterize the ideals of *R, confirming that they are of a simple (non- 
pathological) nature, due to the initial choice of the very well behaved family 
of little-oh polynomials. Finally we will show that, in spite of the presence of 
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infinitesimal numbers h <G *K/o such that h 2 = 0, we can define powers h p for 
every p € R>o and hence, we have arbitrary roots with several good proper- 
ties. This dialectic, between standard structures and the new ring *R, aims 
at presenting the theory of Fermat reals to a general mathematical audience, 
hoping that this will contribute to further studies of this interesting ring with 
infinitesimals. 

2 Metric structures 

The topology used to prove the above mentioned Fermat-Reyes theorem, the 
key theorem for the development of differential calculus of smooth functions 
defined on open sets, is that generated by extensions 'U of open sets U C W 1 . 
In this approach, a subset A C *M" is open if it can be written as 



We will call the resulting topology the Fermat topology. In the present section, 
we want to show that in the ring *R it is possible to define two interesting 
(pseudo) metric structures, corresponding to two different topologies. The first 
one is the Fermat topology, which can roughly be described as the best topology 
for sets having a "sufficient amount of standard points", like, e.g., *U. This 
connection between Fermat topology and standard points can be glimpsed by 
saying that the monad of a standard real r £ K, i.e. 



is the set of all the points which are limits of sequences with respect to the 
Fermat topology (which is not Hausdorff). 

However, in sets of infinitesimals, like the ideal D a , there is only one standard 
point, and indeed the best topology to study this kind of sets is not the Fermat 
one. Therefore, we will define a metric generating a finer topology, called the 
omega topology. When restricted to /i(r), the omega topology is naturally tied 
with the equality up to fc-th order infinitesimals, i.e. x =k y iff °x = °y and 
oj(x-y) < k (see [HIE] for some properties and applications of this notion). It is 
worth noting that the equivalence relation is tied with differential calculus of 
smooth functions defined on infinitesimal sets like D a . Indeed, in [5] it is proved 
that for a smooth function / : D a — > *R there always exist mi, . . . , m n E *R, 
n := [a] being the integer part of a, such that 



A = {J {'U C A | U is open in W 1 } . 



H(t) := {x e *M | °x = r} = {r + h \ h e D^} , 




n 



Moreover, these rrij are unique up to fcj-th order infinitesimals 




where the kj are defined by 
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J 



+ 



= 1. 



o+l 



This permits to define derivatives of smooth functions defined on infinitesimal 
sets. Therefore, it is worth noting that two standard metrics (i.e. with values 
in R and not in *R) are strictly related to the calculus of two different classes 
of smooth functions on the ring of Fermat reals. 

To motivate the definition of our metrics on *R, we can say that: 

• We want to measure the distance between x, y £ *R on the basis of lo(x— y) 
and \°x — °y\ only. 

• We want to extend the classical metric on the reals R. 
Definition 1. Let x, y £ *R, then 



Obviously, both dp and d w extend the usual metric on R; moreover, if x, 
y £ A t ( r )i then d u (x,y) < k iff x =k y. Finally, as we will see later, the idea to 
use also the z-th orders u>i(x — y) to define other metrics on *R is not a successful 
one. 

Theorem 2. The Fermat and the omega metrics verify the following properties: 

1. dp ■ *R x *R — > R is a pseudometric. 

2. i:'lx'I — >R is a metric. 

3. The d u -topology is finer than the dp-topology. 

4- The topology generated by dp is the Fermat topology. 
5. dp and d u are not topologically equivalent. 

Proof: The proof of [7] is direct. Concerning [H we have that d u] {x,y) = iff 
°x = °y and uj(x — y) =0. The order uj(x — y) is zero iff x — y € M, i.e. 
x — y =: r £ R. From °x = °y it follows that r — and hence the conclusion 
x = y. From the definition of order ui{— ), the property ui(— z) = oj(z) follows, 
and hence d^ is symmetric. To prove the triangle inequality, we introduce the 
following lemma, which is a generalization of an analogous result already proved 
for x, y £ Doq only (see e.g. [13] )■ 

Lemma 3. Let x, y £ *R, then 

1. x + y R =>■ lu(x + y) = max [uj{x),w{y)] 



1. d F (x,y) := \°x-°y\ 



2. d u (x,y) := \°x-°y\+u(x-y). 
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2.x + yeR uj(x + y) = i). 

Therefore 

lo(x ~ y) < uj(x — z) + u)(z — y) Vz g *R. 
Using this lemma we have 

d u (x, y) = \°x - °y\ + w(x - y) < 

< \°x - °z\ + \°z - °y\ + lo(x -z) + w(z - y) = 
= d u (x,z) + d UJ (z,y). 

Property [3] follows directly from the inequality d^(x 1 y) < d u] (x,y). 
To prove[7]let us firstly consider a e?F-open set U. Then, for every we can 

find r G R>o such that 

B r (x; d F ) = {y 6 "M | \°x - °y\ <r}CU. (2.1) 

Let U := (°x - r, °x + r) R = {s G M | | x — s\ < r}, then for every y e *U 
we have °y £ U and hence j/ £ W from (|2.ip . Therefore, .t G *U C that is 
W is also open in the Fermat topology. Vice versa, let U be an open set in the 
Fermat topology, then for every igWwe can find an open set U of R such that 
ie'1/C U. Therefore, °x eU and 

(°i - r, °x + r) H C [/ (2.2) 

for some r € R>o- So, for every y £ B r (x; dp) we have \°x — °y\ < r and °y E U 
from (|2.2p . This implies y G 'U CW and proves that B r (x;dF) CM. 

To prove[5|we can consider W = S r (0;d w ), where r > 1. We want to show 
that every ball i? s (0; g?f) is not contained in IA, that is 

3a; g *R : < s , | Q rc| + w(as) > r. 

To prove this it suffices to show that uj{x) > r for some infinitesimal x g -Doo, 
which is trivially true: we can take, e.g., x = dt r whose order is Lo(dt r ) = r 
because r > 1. □ 
Proof of LemmaQH Let Sx := x~°x, Sy := y — °y be the infinitesimal parts of 
x, y, so that, directly from the definition of order, we have ui(x+y) = ui(5x+5y). 
If Sx + Sy — 0, then ui[x + y) = uj(0) = 0. Otherwise, Sx + Sy =/= and from 
Theorem f 2 of |13| we have tu(Sx + Sy) = max [lj(Sx), uj(Sy)] = max [ui(x),w(y)] . 
Finally, it suffices to note that 

Sx + Sy = <^=> a; + y g R 

to prove both properties [7] and [H 

To prove the stated inequality, we note that if x + y g R, then uj(x + y) = 
< oj(x) + u>(y) because w(— ) > 0. If a; + y ^ R, then either uj{x + y) = oj(x) < 
u)(x)+w(y) or uj(x + y) = u>(y) < u>(x)+(d(y). From this, the conclusion follows, 
because lj(x + y) = uj \(x — z) + (z — y)]. □ 
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Definition 4. We will call cu-topology the topology generated by the metric d u . 
It can also be called the topology of the order function (to distinguish it from 
the topology of the order relation) . 

Is it possible to generalize using higher orders Wj? 

It is very natural to try a generalization of the metric d w considering, e.g., also 
the information given by to 2 (x — y): 

d 2 (x, y) ■= \°x - °y\ + oj(x -y)+u: 2 (x- y). 

However, an immediate problem is that higher orders u>i(— ), i > 1, do not share 
the good properties of w(— ). For example: 

• If x = dt and y = dt^ + dt 2 , then u) 2 (x + y) = ui 2 {y), but if y = dt 2l then 
uj 2 {x + y) = lo(x). 

• If u>(x) = u(y) and °x\ + °y 1 =/= 0, then co 2 (x + y) = max [u)2(x), w 2 {y)\, 
but if °x\ + °y 1 = 0, uj 2 (x) = uj 2 (y) and °x 2 + °y 2 ^ 0, then uj 2 (x + y) = 
m.<nt[wz(x),U3(y)]. 

• If x = d^6 + 2 dt , y = dt-? + dt and z = 3 dt , then uj 2 (x — z) + uj 2 (z — y) < 
u 2 {x — y). Therefore, the inequality of Lemma [3] cannot be proved for 
higher orders. 

We will solve this problem with the following result. 
Proposition 5. Let i G N>i, and suppose that 

i 

di(x,y) := \°x-°y\ +^a>j-(x- y) 

3=1 

verifies the triangle inequality. Then di and d u are equivalent. 

Remark 6. In this statement, we mean U3j(r) := for r € R and 0Jj(x) := if 
j > N, where N is the number of summands in the decomposition of x. 

Proof: It is readily verified that the di-topology is finer than the w-topology. 
To prove the converse, let us first consider y £ B r (x;d ul ), with r < 1. Then 
du(x,y) = \°x — °y\ + uj(x — y) < r < 1. Therefore, uj(x — y) < 1 and hence 
uj(x — y) = 0. This means i-i/el and hence u>j(x — y) = and d,(x, y) = 
d LJ (x,y) < r. Therefore, B r {x]d UJ ) C B r (x;di) if r < 1. Finally, for a generic 
ball B s (x;di), take n £ N>o such that — < 1, then B s (x;di) D B s / n (x;di) D 

Bg/n^du). □ 

2.1 Fermat and (^-completeness of *R 

The proof of the following result follows directly from the equality d~F{s n , s m ) = 
I $n I ■ 
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Theorem 7. With respect to the Fermat metric dp, the ring *R is complete. 
In particular, if (s n ) n is a Cauchy sequence with respect to dp, then (°s n ) n is a 
standard Cauchy sequence of R. Let r G R be its limit, then 

V.x G : {s n )n converges to x with respect to dp. 

Before studying the w-topology, we want to understand better the intuition 
underlying this metric, because it is strictly related to nilpotency of every in- 
finitesimal of *R. Let us start with an example: 

s n := — h dtjv Vn G N>o, 
n 

where N G N >0 . We have d u (s n , 0) = £ + u (jj + dtN - o) = ± + f ->• as 

n — > +oo. Therefore, (s n ) n converges to in the w-topology. However, exactly 
because of nilpotency, we also have 

Vn > N : dijv = hence s n = — G R. 

n 

More generally, if d u (s n ,x) — > 0, then |°s„ — °x| — > 0, but also w(s„ — x) — ?> 
0. This means that the order ui(s n — x) goes through smaller and smaller 
infinitesimals. However, because of nilpotency, infinitesimals of "R cannot have 
order less than 1, and hence the order ui(s n — x) must collapse from 1 to 0. The 
following theorems will clarify this intuition. 

Theorem 8. Let x G *R and s G (0, 1]r ; £/ien 

-B s (z; d u ) = (s, s) K + {a;}, 

i/iai is, every y G d w ) can fee written as y — r + x, with r G R, — s < r < s. 

Proof: For y G B 8 {x; d w ) we have \°y — °x\ < s and ui(y — x) < s < 1, therefore 
— x) = and j/-i=:rgl. Moreover, r = °y — °x so that |r| < s, that is 
r G (— s, s)r. Vice versa, if y G (— s, s)r + {x}, then y = r + x, with |r| < s. So 
|°y — °x\ = \r\ < s and uj(y — x) = uj{r) = and d u (y, x) < s. □ 
The following theorem characterizes w-convergent sequences formalizing the 
intuition presented above. 

Theorem 9. Let (s„)„ fee a sequence o/*R, then we have that 

lim s„ = x G *R (2.3) 
mt/i respect to the omega topology if and only if the following conditions hold 

1. lim„^ +00 °s n — °x in R 

2. The sequence of infinitesimal parts is eventually constant and equal to Sx, 
i.e. 

3N G N : Vn > N : Ss n = Sx. 
(where we recall that Sy := y — °y.) 
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Proof: We only have to prove [H as the rest of the proof is immediate. Because 
of (|2.3[) , we have that 

lim u>(s n — x) = = lim uj (Ss n — Sx) . 

n— >+oo n— >+oo 

Therefore to (Ss n — Sx) < 1 for n > N, and hence u (Ss n — Sx) = 0. This means 
that Ss n — Sx € R and hence Ss n — Sx = because Ss n , Sx £ Doq. □ 

Example 10. s n = — + dtw + dt-2 — > dt2, whereas s„ = dt 1+ i_ is not an 
w- convergent sequence. 

Theorem 11. The ring *R is complete with respect to the metric d^. 
Proof: Let (s n ) n be an w-Cauchy sequence, then 

lim \°s„ - °s m \ = 0, 

n— >+oo 
tn— v+oo 

so that the lim„^ +oc °s n =: r £ M exists. Moreover, we have 
lim u>(s m — s n ) = = lim cu(Ss m — Ss n ). 

n— >+oo n— >+oo 
m— > + oo m— )- + oo 

Like in the previous proof, we have 5s n = 8s m = 8s n for n, m> N. Therefore, 
setting x := r + 5sjv £ *K we have lim n _j. +00 s„ = x. □ 

2.1.1 The omega metric on Doo can be defined by a pseudovaluation 

In this section we want to prove that the restriction of to the set of all the 
infinitesimals 

duih, k) = uj(h — k) V7i, k s £*oo 

is induced by a pseudovaluation. We will see that the same idea doesn't work 
outside Doo . These notes have also the aim to fix some small error made on the 
same topic in [T3] , 

For completeness, we start from the definition of pseudovaluation on a 
generic ring with values in RU {+oo}. 

Definition 12. Let A be a ring, then we say that v : A — > M. U {+00} is a 
pseudovaluation if the following conditions hold: 

1. v. A — > R U {+00} 

2. Vx e A : v(x) = +00 <=^> x = 

3. v(x ■ y) > v(x) + v{y) 

4. v(x + y) > min [v(x),v(y)] 

5. v(x — y) = v(y — x) 



10 



Obviously, the last three conditions are supposed to hold for every x, y G A. 
Remark 13. 

1. We assume the usual rules about the relationships between +00 and order 
or sum on R: +00 + r = +00 > r for every r € R. 

2. The essential difference between a pseudovaluation and a valuation is prop- 
erty [3J where an inequality replaces an equality. Indeed, it is not hard to 
prove that equality cannot hold in any ring with zero divisors, like *K. 

To motivate the necessity of our definition of pseudovaluation, we anticipate 
the following lemma, with which we can treat terms of the form w{x ■ y). 

Lemma 14. Let x, y S *R, then 

lo(x) if x G Doo, y g Drx, 

^k + ^kY 1 tf*,veA»\{0} 

Moreover, if x, y ^ R, then 

u(x ■ y) < uj(x) ■ u(y). 

Proof: Using the infinitesimal parts, we can write x ■ y = °x ■ °y + °x ■ Sy + 
°y ■ Sx + Sx ■ Sy, so that, by the definition of order, the equality ui(x ■ y) — 
lu(°x ■ Sy + °y ■ Sx + Sx ■ Sy) follows. 

In the last case of the statement, i.e. x, y £ Z?oo \ {0}, we have 



uj(x ■ y) 



lu(x ■ y) = uj(Sx ■ Sy) = 



1 



uj(x) u(y) 



u(x) ■ uj(y) 
u(x) +ou(y) 



< 



< -u)(x) -u(y) < u)(x) -uj(y) 



Moreover, let us also note that 



j(x-y) 



1 

o(x) 



and hence 



u(x-y) <u{x). (2.4) 

Below we will use this inequality for suitable infinitesimals x and y. 
li x € Doo and y £ -Dod, the cases Sx = or Sx ■ Sy — are immediate. 
Otherwise, we have 

oj(x ■ y) = ui(°x ■ Sy + °y ■ Sx + Sx ■ Sy) = to{°y ■ Sx + Sx ■ Sy) = 

= max [oj(°y ■ Sx),ui(5x ■ Sy)] . (2-5) 

Now, we can apply (|2.4p to the product Sx ■ Sy, obtaining 



ui(Sx ■ Sy) < tu(Sx) = uj(°y ■ Sx). (2-6) 

Note that the last equality is due to the fact that y ^ Doc, so that °y ^ 0. 
Therefore, from (|2.5p and (|2.6[) we obtain oj(x ■ y) = ui(°y ■ Sx) = uj(Sx) = oj(x). 
Finally, ui(x) < u(x) ■ to(y) if y ^ R because in that case oo(y) > 1. □ 
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Remark 15. In the statement of the previous theorem, the case where x, y $ 

is not included. This is done to avoid an overcomplicated statement. Indeed, 

we have several sub-cases: 

1. If °x ■ 6y + °y ■ Sx + Sx ■ Sy = =s> u(x ■ y) = 



2. If °x ■ Sy + °y ■ Sx = and Sx ■ Sy ^ =^> uj{x ■ y) = 

3. If °x ■ Sy + °y ■ Sx =/= =>■ w(a; ■ y) = max w(y)] 



i i i 



The idea to define a pseudovaluation is to derive the property v(x • y) > 
v(x) + v(y) from u(x ■ y) < ■ w(y), so that the natural try is the following 

Definition 16. For every infinitesimal h € £>oo, we define 

v(h) :=-]og[u{h)], 
where we use the convention that — log(O) := +oo. 
Remark 17. 

1. The metric associated to v is e ~ vlyh ^ k ^ = uj(h — k) = d u {h, k). 

2. If we define v(x) := — log + lo(x)] for every x £ *R, we do not obtain 
a pseudovaluation. Indeed, we have e.g. w(3 + 5) < min [w(3), w(5)]. 

Theorem 18. v : Doo — > M U {+00} is a pseudovaluation on the subring 
(ideal) Doc of all the infinitesimals. 

Proof: We only have to prove property 0] of Definition [T2J as the others are 
immediate. We can suppose x + y 7^ 0, because otherwise the proof is obvious. 
Since x, y £ and i + we have that lu(x + y) = max [u}(x),lu(jj)]. We 

will proceed in the case uj{x) > uj{y), the opposite being analogous. Therefore, 
oj(x + y) = uj(x) and v(x + y) = — log [w{x + y)] = v{x). >From w(x) > uj(y) we 
have — logw(a;) < — logw(y), so that v{x) < v(y) and min [v(x), v(y)] = v{x) = 
v(x + y), which is our conclusion. □ 



3 Ideals and their characterization 

In this section, we want to study the ring of Fermat reals from the point of view 
of some standard algebraic structure. To begin with, we note a few elementary 
algebraic properties of *R: 

• There are no nontrivial idempotents in *R. 

• *M is not reduced. 

• x £ *R is a zero divisor iff a; £ -Doo iff x is non-invertible. 
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• *R is an exchange ring (i.e., for each x £ *R there exists an idempotent 
such that x + e is invertible) . 

• *R is an /-ring (lattice ordered ring), as well as a normal /-ring (cf. [2]). 

Next we want to study the ideals of *R. We will see that, as opposed to 
other rings containing infinitesimals (see, e.g., |27] ) the ideals in *R can be ex- 
haustively described. Of course, this is essentially due to the very simple family 
of little-oh polynomials used as representatives of new numbers in *R. We start 
by proving that the only maximal ideal is the set of all the infinitesimals. 
The idea to consider only "well behaved" functions x : R>o — > R (the little-oh 
polynomials), in the definition of the ring *R, is tied with the fact that Fermat 
reals do not represent a new foundation for the entire calculus. Indeed, our aim 
is only to extend ordinary smooth functions, so that it suffices to evaluate them 
on 'well-behaved numbers'. 

The situation is entirely different in NSA, which aims to be a new, indepen- 
dent foundation of the whole calculus. For example, suppose we want to prove 
that ordinary continuity of a function / : R — > R at xq £ R is equivalent to 

Vx£*R: x~x =>■ *f(x)~*f(x ), (3.1) 

where x ~ y means that x — y is infinitesimal. However, (|3.1|) is nothing more 
than the continuity of the function / stated using sequences, i.e. 

VieK N : 3 lim x n = x =^3 lim f(x n ) = f(x ). (3.2) 

This is equivalent to ordinary continuity only if it is stated for every sequence 
x £ R N , as is obvious from the corresponding proof. 

Lemma 19. Let I be a proper ideal of the ring *R ; then I C Dqo. 

Proof: Let x £ I and suppose, by contradiction, that °x ^ 0, then x would be 
invertible and for every y £ *R we could write y — y ■ x~ x ■ x. By hypothesis 
x £ I, which is an ideal, so we would have y £ I, that is I = *R which is 
impossible because / is proper by hypothesis. □ 
Directly from the decomposition of every Fermat real it follows that 

•R/Doo ~ R 

and hence is a maximal ideal. With the following result, we prove that in 
fact it is the only one. 

Theorem 20. Let L be a proper maximal ideal of the ring *R ; then L = D x . 

Proof: By the previous lemma, we have that / C D^,. The ring *R is commu- 
tative and with unity, so *R// is a field because, by hypothesis, I is maximal. 
Take x £ D^. We have two cases: either x + I = I or x + Iy^L. In the first 
one, we have x + = x £ I. In the second one, x + 1 would be invertible in the 
field *R/I, so, for some y £ *R, we can write 

(y + I)-(x + I) = l + I i.e. yx + I =1 + 1, 
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that is yx + i = 1+j for some i, j £ I. Taking the standard parts in this equality 
we obtain °y ■ °x + °i = 1 + °j. >From / C we deduce that °i = °j = 
and, therefore, that °y ■ °x = 1, that is °x ^ 0, which is impossible because 

ieD m . □ 

This proof can be easily generalized to the following 

Theorem 21. Let R, F be two commutative rings with unity, with F non 
trivial. Moreover, let s : R — > F be a ring morphism and set D := ker(s). 
Finally, let us suppose that 

\/x £ R : x D => x is invertible in R. 

Then D is the only maximal ideal of R. 

In our case s = °(— ) is the standard part map. 

Finally, we prove that every proper ideal is either of the form 

D a = {x G Doo | oj(x) < a + 1} for some a € R>x U {oo}, 

or of the form 

I a := {x 6 Doo | < a} for some a £ M>i. 

As we mentioned above, the first type of ideal is used in infinitesimal Taylor 
formulas, whereas the second is used in the study of the equivalence relation = Q 
of equality up to a-th order infinitesimals. 

This characterization is tied to the possibility to solve in 'M. the following 
class of linear equations. 

Theorem 22. // a, b, c £ *K and a < c < a + b, then 

3x £ "R : a + x - b = c 

For the proof of this theorem, see |101 [5]. Let us note that we cannot have 
uniqueness of solutions, due to nilpotency. For example, if a = 0, c = d^2 + dt 
and b = (H3, then x = dt§ + d£ 3 / 2 is a solution of a + x ■ b = c, but x + dt is 
another solution. Moreover, let us note that this theorem is not in contradiction 
with the non Archimedean property of 'M. (let a = and b £ Doo) because of 
the inequalities that c must verify for a solution to exist. 

Using this result, we can prove the desired characterization: 

Theorem 23. Let J be a proper ideal of *M, and & £ J. Moreover, set 
O(J) := £ M> | j £ J} , a:=supO(J). 

Then 

1. Vc £ 'R : -b < c < b => c£j 

2. a £ 0{J) =*> J = I a 
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3. a<£0{J) a>l and J = D a _ x . 



Proof: To prove [71 set a = in Theorem |2"2"1 Then since — b < c < b, we can 
distinguish two cases (we recall that the order relation in *R is total). If c > 0, 
then < c < b and we can hence solve the equation x ■ b = c and, therefore, 
c G J because b E J . Otherwise, c < and so < — c < 6. We solve the 
equation x ■ b = — c, that is (— ic) ■ b — c so that c € J again. 

To prove we first note that J C I a by the definition of a. Vice versa, let 
i E I a , i.e. < a. By hypothesis, a £ O( J), so that we can write a = cu(j) 
for some j <G J. We can suppose j > because, otherwise, < -j € J and 
a = uj(j) = uj(—j). We distinguish two cases. If oj(i) < a = cj(j), then i < j 
by the properties of the order relation we mentioned in the introduction (see 
Theorem 4.2.6 in jH]). On the other hand, we also have that w(— i) = u>(i) < uj(j) 
and hence — i < j because j > 0. Therefore, —j < i < j, and the conclusion 
i <E J follows from[7J Let us note that, in general, we have just proved that 

VijG'R: j>0, u{i)<u>(j) -j<i<j. (3.3) 

In the second case, we suppose that oj(i) = a so that, by the decompositions of 
i, j and for suitable a, ft £ R^o and h, k € -Doo, we can write 

j = a ■ dt a + h , uj(h) < a (3.4) 
i = ft ■ dt a + k , ui(k) < a. (3.5) 

Therefore, from p.4p and Q3.3P it follows that —j<h< j and hence h E J 
from property [JJ So a ■ dt a = j — h € J and df Q € J because a^O. Hence, 
we also have that ft ■ dt a e J. Finally, cj(fc) < a = oj(j) so that —j<k< j 
from (|3.3p and k G J from property [T] We have proved that ft ■ dt a , k E J, so 
i = ft ■ dt a + k G J, which is our conclusion. 
Finally, to prove [3] we first note that 

D a -i = {x E Doo | lo(x) < a} , 

where we use the conventions oo ± 1 = oo. If a ^ O(J) then we have < a 
for every j G J, and therefore J C D a _i, considering also that J C Z?^. 

Vice versa, if i E D a _i, then we can find j E J such that uj(i) < ui(j) < a 
because a = sup O(J). As usual, we can suppose j > 0. From (|3.3p . it follows 
— j < i < j and hence i E J by property [JJ To finish, let us note that because 
^(j) > 1 or Lo{j) = we necessarily have that a = or a > 1. However, the 
first possibility would imply J = {0} and hence a G O(J), which is impossible 
by hypothesis. □ 

4 Roots of infinitesimals 

In the ring of Fermat reals *R, the existence of non zero nilsquare elements: 
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h ^ , h 2 = 0, 



(4.1) 



is incompatible with the existence of a square root and of an absolute value with 
the usual properties. In other words, if we want to define roots of infinitesimals, 
we have to avoid from (|4.ip the following inference: 

h 2 = therefore Vh? = Vo = 
Vh 2 = \h\ = therefore h = 0. 

We recall that only smooth functions / : R — > R can be extended to *R. 
In particular: 

• Because they are locally Lipschitz, these functions verify 

Vx, y G *R : x = y in *R =s> fox = foy in *R. 

• Because they are smooth, they take little-oh polynomials into little-oh 
polynomials: 

V.t G K„[t] : /ox G R Q [i]. 

It is hence natural to expect some problems extending, e.g., the square root to 
the whole of *R. 

The first natural solution is to extend the roots only where they are smooth, 
i.e. on R^o- This is equivalent to defining the roots only for invertible Fermat 
reals (and positive in case of even roots or irrational powers) . For details about 
this approach, see [5], section 4.3, or [TT], section 12. 

Another problem we have to take into account, and concerning roots of 
infinitesimals, is that the equation x 2 = c, for c G -Doo, always has infinitely 
many solutions, e.g. 

(dt 4 ) 2 = dt 4 / 2 = dt 2 

( dt 4 + hf = dt 2 + h 2 + 2hdt 4 = dt 2 V/i G : u(h) < -. 

3 

Therefore, we have infinitely many square roots of an infinitesimal. This 
means that, although in R we have that (— ) 2 : R>o — > R>o is bijective and 
: R>o — > R>o is its (left and right) inverse, in *R>o we don't have injectiv- 
ity, and, therefore, we can have, at most, a right inverse. Indeed, we will prove 
that (— ) 2 : *R>o — > *R>o is surjective, and 

Vx G *R>o : (a/^) 2 = x but 3k : \fx 1 ~ =k x. 

Because we have infinitely many solutions of equations of the type x p = c, 
p G R>i, a first idea is to choose, among them, the simplest solution. Here, with 
"simplest", we mean "the solution x without unnecessary terms in its decompo- 
sition, i.e. without terms that become zero taking the power x p, \ A similar 
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idea of "simplest solution" has already been used in [§] to define derivatives of 
smooth functions defined on infinitesimal sets. For example, both x = d<4 and 
x = dt± + dt are solutions of the equation x 2 = dt2, but, intuitively, the first 
one is simpler compared to the second one, which contains the unnecessary term 
dt. 

However, there is another, more manageable idea to define roots of infinites- 
imal numbers. Let 

N c 
i=l 

be the decomposition of c £ £>oo- Suppose c > 0, so that lu(c) = lo\{c) > u>2(c) > 
• • • > ujpf c (c) > 1 and °Ci ^ 0, c\ > 0, then, for p £ M>o, we would like to write 





%=2 C1 / 



■dt 



(4.2) 



where 



oo 



a © b 



aQ b := 



Va,6e 



However, the right hand side of (|4.2j) is now a well defined term, because the 
base of the p-th power is invertible. 

Remark 24. 

1. Note that the right hand side of (|4.2[) is well defined if °c\ ^ 0, i.e. if 
c £ Doo \ {0}, and because ui(c) > w,(c), so that U)i(c) G w(c) > and 
hence d£ Wi ( c ) ew ( c ) is well defined. Moreover, it is not hard to prove that 
LUi(c) u(c) > 1 if i > 1. 

2. It can be useful to note that setting 

Qb := -b 
a © b := a ■ b, 



we easily have that (Roo, ffi, G, 0, oo) is a ring and the reciprocal function 
jpr : Rqo — > R is a ring isomorphism. 
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Definition 25. Let c <E £>oo, c > 0, and p € K>oj then 



:= (° Cl ) p dt^-il + Y^^ dt m(c)eu(c) \ . 

Of course, if p = ^, where m, n £ N and n is odd, the hypothesis c > can be 
dropped. 

Example 26. 

1 . Let us find \f~d~i using the previous definition. In this case, we have N c = 1 , 
°ci = 1, uj(c) = 1, so 



dt=(l) 1 / 2 dt 1/2 -(l + 0) 1 / 2 = dt 2 . 



2. We want to find V dt 2 + dt: 



V dt 2 + dt = dt 4 • (1 + d<i e2 ) 1/2 = dt 4 (1 + dt 2 ) 



We recall that 



>l/2 



dt 4 • ( 1 + 



2 /is 

2 1 dt 2/r 

n 1 



n=l 

1 



1 



dt 4 ■ ( 1 + - • dt 2 - - ■ dt 



dt 4 + -dt 4/3 . 



(1 + /if = 1 + ( • k fl VheD 



n=l 



(4.3) 



because of the elementary transfer theorem (Theorem 24 in Q]T|). Finally, 
let us note that the series in (|4.3p is really a finite sum because of nilpo- 
tency of every infinitesimal h £ Z?oo. 
As expected, we have that ( dt 4 + ^ dt 4 /^ = dt 2 + dt. 



3. yj dt 2 - dt 3 / 2 - dt = dt 4 - \ dt 12 / 5 - | dt 12 / 7 - ^ di 4/3 - dt 12 /u- 
Generalizing these examples, we have that 

+oo 



'° cl f dt. 



1 + H ( ! ) ' ( £ ^77 d ^.( C )e^(c) 



In the following theorems, in considering x p for generic p G K>o, we will 
always suppose x > 0. However, this hypothesis can be dropped in case of odd 
roots, and the proofs will remain essentially the same. 
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Theorem 27. Let x £ -Doc; x > 0, and p £ R, < p < 1, t/ien we /iawe: 
-Z. (a; p ) p = x 

2. If X? 7^ and k := max ^^(x), u>2 }> ^ erl ( a; '') ~ fc ■ r - 

Remark 28. 

1. To understand better, it can be useful to clarify what is the difference 
between the computation of (x p ) 1 ^ p and that of (x 1 / p ) P : 

(a) (x p ) 1 ^ p : Because p £ [0, 1]r, the computation of x p is included in the 
Definition 1251 Therefore, we must: 

i. Express x using its decomposition. 

ii. Use Definition l25l 

iii. With the obtained result, we finally have to compute the subse- 
quent power (— However, | > 1, so that this operation is 
smooth and doesn't present any problem. 

(b) (x 1 / p ) P : In this case, the situation is the opposite one. 

i. jj > 1, so the operation x x l p =: y is smooth. 

ii. However, to compute y p , we must apply Definition I25[ so we 
firstly need the decomposition of y = x x l p . Of course, it is not 
easy to find this decomposition as a manageable function of the 
decomposition of x. 

2. Let us note that if x, y £ Drx, and k := max [oj(x), oj{y)], then it is trivially 
true that x y, because ui(x—y) = k if x ^ y, and w{x~y) = otherwise. 
This shows that [H of Theorem [27] is not trivial. 

3. Theorem [27] represents an overcoming of the incompatibility between nil- 
potent infinitesimals and existence of roots. Indeed, property can be 
applied only if x^ p ^ 0. Moreover, if h £ D x \ {0} and h 2 = 0, we have 



h 2 = V0 = 0, but, in general, Vh 2 ^ \h\, e.g. \J ( dtf = ^ | dt| = At. 

Proof of Theorem 12 Tt Let x = J2f=i bj dtp j be the decomposition of x. 
Because < p < 1, from Definition [25] we have 

/ N ^ 

x p = b p dt ix ■ 1 + V/dt fte/3l 
P V 3=2 bl 

Now, we have to apply the power (— ) l ' p , which is smooth and has the usual 
properties of powers (see, e.g., [SJ, section 4.3). Therefore, we can write 
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(x p )v = Kdta 



AT 



J=2 



j=2 



This proves [I] 

To prove we firstly have to compute the smooth power (— ) 1 / p 

JV 

The idea is to use the usual properties of (— ) 1 ' p and to gather up the leading 
term b\ dtp 1 : 



x '• 



N 



i=2 

AT 



3=2 



(4.4) 



We are not able to find the decomposition of this number, but we can surely 
claim that 



— p ■ Pi — p ■ lj(x) if xv ^ i.e. if p ■ Pi 



> 1. 



(4.5) 
(4.6) 



This guess is based on the idea that in (|4.4[) . the infinitesimal 6^ • dip^ is 
multiplied by an invertible number, whose standard part is 1. Indeed, we have 
the following 

Lemma 29. Let k, h G and y £ *R such that y is invertible and k = h ■ y. 
Then 



uj{k) = u>(h) 
°ki = °hi-°y. 



20 



We postpone the proof of this lemma to the end of the current proof. 
Using (|4.5[) and (|4.6[) and Definition [25l we can write 

(M ( X p\ \ P 

l + Y.^r J±dt ( x \ , A . ( 4J ) 

where M is the number of terms in the decomposition of x x l p (of which, we 
really know only the first term) . Using again Lemma [29] applied to ()4.7j) , we 
have 



U>(x). 



Let us observe that x > 0, hence x ^ and ui(x) > 1, so that ui (°xx dt^M*) = 
uj(x). We have hence proved that the first terms in the decompositions of both 

x and ( X? } are the same. Therefore 



< max 



0J2 



^xp^j | = k 



and hence . □ 

Proof of Lemma I29t Because y is invertible, we have that °y ^ 0. Write the 
product k = h ■ y using decompositions 



N k / N h \ I N y 

k = Y^ ° k * d ^,(fc) = I X ° h v dt u P (h) • \°y + X) % d W'«) 

= °hi°ydt u ( h ) + ^2°h p °ydt UJp{h) +^^°h p °y q dt u]p[h)mulq(v) . (4.8) 

p=2 P=l <3=1 

However, in general, we have a © & < min(a, 6) for every a, b G R>o, so 
that, in (|4.8[) the leading term is °h\°y dt w ^ and hence from the uniqueness of 
decomposition, the conclusion follows. □ 

Remark 30. 

1. Let us observe that, in the hypothesis of Theorem [27j we also have 

y = x p =*> (yp^j =y. 

ii / i\p 

In fact, yp = {x p ) p = x, and hence [y p ) = x p = y. Therefore, for 

numbers of the form y = x p , the equality fj| of Theorem 1271 becomes exact. 



21 



We can interpret this result saying that our Defmition[25]of c p gives exactly 

i 

the simplest solution of the equation x? = c. Indeed, like in the case 



(d<2 + dt) =1 dt 2 + dt, we can say that in[J|we don't have an exact 
equality if the number x contains unnecessary infinitesimals with respect 
to the power (—) 1 ^ p , like dt in the previous example. See section [4731 for a 
formalization of the notion of "unnecessary term with respect to the power 
(_)i/p» 

2. The equality [H up to infinitesimals, implies that 

*R>o — > *R>o is not surjective. 

For example y = dt2 + dt cannot be written as the square root of some 
number x. Otherwise, we would have 

y = y/x V 2 = ( dt 2 + dtf = dt = (V^) 2 = x, 

but then *Jx = ydi = dt 2 ^ y. Of course, this corresponds to saying 
that the square is not the right inverse of the square root. 

3. Trivially, we can consider a smooth function g : R — > R having a root of 
order n G N>o at x = 0, i.e. such that 

g{h) = a-h n Mhe D n , 

where a <E R^o- We can hence define a sort of infinitesimal right inverse 
of g, setting 

f{h) := y| Vfee 
if n is odd or n is even and a > 0, and 

f(h) := { — VhtDoc 
V —a 

if n is even and a < 0. Then we have g(f(h)) = ±h for every /i G -Doc, 
with the positive sign in the first case. 

4.1 A formula to compute a root 

By definition, if c G £>oo \ {0}, we have 

:= (° Cl f di^ • ( 1 + ^ ^ d Wc)ew(c) . 

The p-th power of the invertible term can be computed in several, obviously 
equivalent, ways. 
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1. Using the infinitesimal Taylor formula of the function (1 + x) p , with 



i=2 



Cl 



(4.9) 



2. Equivalently, we can use the formula (1 + x) p = 1 + Yln=°i (n) ' * T "> ^ or 
| a; | < 1, which transfers to by Theorem 24 of |13] , 

Applying the second method, we get 



C p = (° Cl ) p dt^ 



(°d) p dt^ 



71=1 



£ ST" d ^i(c)ew(c) 

v 1=2 1 , 



n\ Nc ( °Ci x 7<_I 

n ( o„ 

|7|=n 



dt^(o 



where 



if p £ N 



min ( r cj 2 (c) — aj(c) n ,p) if p G N 



(4.10) 



with r a n G N the ceiling of a G R, that is the smallest integer greater than or 
equal to a. Note that the first alternative of (|4.10p is due to (|4.9p . whereas the 
second one is also due to the equality = if n > p. 
Using Theorem 13 of [13], we have 



T dt „ iMe „ M ? o V 

1-i — 1 

i=2 i=2 



7i-i 



w<(c) w(c) 
Wi(c) G w(c) 



< 1 



> 1. 



7i-i 



In the following, we will set w(c, 7) := (B^ 2 ^'^^ 9 ^ c ^ , so that 



l= 2 



(c,7)" 



Finally, we obtain the formula 



fe c . 



„7i 



>n/ z — ' 7! 

n=l v 7 7er c ,„ ' 



o_7-W c -i 1, 



(4.11) 
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where 

r c ,„:= {7 e N^ 1 I | 7 | = n , w( Cj7 )>1}. 

4.2 Properties of roots: the general theorem 

For generic x, y <G *R we can state the following 

Theorem 31. Let x, y € be strictly positive infinitesimals, and p, q € R>o, 
then: 

1. uj [{x p ) q \ = u(x pq ) =: 01 and [(x p ) q } 1 = (x pq ) 1 . 

2. lo [{x ■ yf] = lo{x p ■ yP) =: o 2 and [(x ■ y)\ = (x p ■ y p ) v 

3. 3k e R : 1 < k < 01 and {x p ) q = k x pq . 

4. 3k e R : 1 < k < o 2 and (x ■ y) p = k x p ■ y p . 

5. x p -y q = x p+q . 

Before proving this theorem, we need the following very useful lemma: 

Lemma 32. Let c = ^2jLi a j °-^ctj > with a\ > ctj > 1 for every j = 2, . . . , M , 
and a\ > 0. Let us note explicitly that not necessarily this is the decomposition 
of c. Then 

c p = a p dt^ ■ h+^^dt^emj • 

Of course, this lemma states that the formula used for the definition of c p 
can also be used starting from a representation c = YljLi a j ^t aj which is not 
necessarily the decomposition of c. To apply this lemma, the important step is 
to find the greatest infinitesimal a\ dt ai and to check that all the other terms 
aj dt aj can be zero only if aj = 0. 

Proof of Lemma I32t Starting from c = 5Zj=i a j , we firstly sum all the 
coefficients aj having the same infinitesimal dt aj , i.e. if 

a q ■= \j = l,-..,M , a, = q} Vg 6 {oj\j = 1, . . . ,M} 

0:={ aj \j = l,...,M , a aj ^0}=:{q u ...,q N }, (4.12) 

then 

N 

c = ^2 a q dt q = ^a 9l dt qt . (4-13) 
qeO »=x 

Let us note that in (|4.12[) , q±, . . . , q^ is any enumeration of the elements of 
the set of all orders O. Now, all the summands in (|4.13|) are non zero, because 
of our definition of the set O. Therefore, reordering the summands in (|4.13p . 
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we obtain the decomposition of c. Formally, this means that we can find a 
permutation a of {!,..., N} such that 



JV 



(4.14) 



is the decomposition of c. Let us note that, to obtain (|4.14p . we need that 
for every i = 1 , ... N we can find j = 1 , . . . M such that qi = aj > 1 . By 
definition of decomposition, q ai is the maximum order in (|4.14p . i.e. q ai = 
max {qi, . . . , qw} = max{aj \aj € O}. However, we have that q ai = a\, be- 
cause «i > ay by hypothesis, and because 

a ai = y~]{a,j\j = l,...,M, aj = a{\ = ai ^ 0, 

so that ot\ G O and so a qtJi = a ai = a\. We can now apply our Definition 1251 
using the decomposition (|4.14[) : 

/ N - \P 

c? = a p 1 dt^ ■ h+^-^ldt^eaxJ • 
Now, we only have to retrace the previous steps, so as to eliminate a, q, a, etc. 



v 



i=2 



( V 

/ 



q eo 

( M 



"1 



M 



1+ ^ — d^Q-ecu + — dt a .Q ai 



J=2 



a i 



3=2 

O a .=0 



a? dt^i • 1 + V^d^e^ , 



which is our conclusion. 



□ 



Proof of Theorem 13 It To proveQJ let x = Ylj=i bj dtp j be the decomposition 
of x. The idea is to use formula (|4. 1 1|) to compute x p , and then Lemma [32] to 
compute (x p ) q . To avoid heavy notations, we will use the simplified symbols 
k ■= k X p and T \— T^ x n i 
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(x p ) q = \b p 1 dt Ex + 

n=i v 7 7er 

Before using Lemma l32j we need to prove that d£^i_ is the greatest infinites- 
imal, so let us compute 




So, we need to prove that -Jr- < + 53i=2 ^at ^ s n < 7»-i ' fr- 

In fact, /?i > /?i so that J^LaT*-! ' f 1 > £^2 7*-! = n - Moreover, we can 
suppose to restrict the set T to those 7 such that oj(x, 7) — > 1, because, 
otherwise, the corresponding term di,^ 7 w£i = 0- We can ?rence apply the 

Lemma [521 obtaining 



{xV) q = b\ q dt^ 



7er 



(4.15) 



On the other hand, we have 



N 



6f di a • l + V-^di^ 



3=2 



(4.16) 



Therefore, the conclusion follows from (|4.15[) . (|4.16[) and Lemma I2U1 

To proved! we can use the same method as before. Let y = X)fc=i e kdt 
be the decomposition of y, then 



JV M 

(xyf = ( ^^6 iefe dt ft( 
, j=i k=i 
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Of course, j3\ e± > j3j for any j and k, and we can also consider 

I := {{j, k) \j = 1, . . . ,N , k = 1, . . . ,M , ft © e k > 1} , 
so that to the sum 

\U,k)ei 

we can apply Lemma 1321 We obtain 

( 



{xy) p = b\e\ dt^ieei 



On the other hand, we have 



1+ X jLL£*d 



V 



(i,fc)e/ 
Cj-,fe)#(i,i) 



6i • ei 



ft-ee fc e(/3ieei) 



= b p dt Sl • I 1 + ^-^d^e^ 



6i 

J=2 




e? dt^ ■ 1 + V ^ di efcQei = &?e? di^ • (1 + hf , 



where h G is obtained from the product of the previous p-th powers with 
invertible bases. Once again, the conclusion follows from Lemma [2U1 

The proofs of [3] and [7] are straightforward, taking into account [TJ and [H so 
that, e.g., in the difference (x p ) q — x pq there appear only infinitesimals of order 
greater than o\ . 

Finally, to prove [5| we only have to apply our Definition 1251 of power: 



N u 

OA 



N \p+i 
i=2 



which is the conclusion. □ 



27 



Remark 33. For generic x, y G -Doo, properties [3] and [7] of Theorem I3T1 cannot 
be improved. Indeed, if we had always 

(x p ) q = x pq Vx G £»oo Vp, g > 

we would have, as a consequence, the general validity of 

(x?y = x? p = x Vx G £oo Vp G (0, 1] R , 

but we know that this property is not generally true. 
Analogously, from the general validity of 

(xy) p = x p -y p V.t e D M Vp > 0, 

we would have 

(x 2 ) 5 = (x ■ i)5 = • = (\/a;) 2 = x Vie Dqo, 

but we know that this is not generally true. 

Therefore, on the one hand these seem the best results attainable. However, 
it does not seem desirable to work with the equality up to infinitesimals of 
some order, in particular for such basic operations. 

We have already noted that several counterexamples are of the form 

i 

(dt 2 + dt) 2 

where we have terms like dt which are, intuitively, unnecessary with respect to 
the square. Our next aim is to formalize the idea of incomplete term with respect 
to (— )p, and to prove that the usual properties of the powers hold, with the 
usual equality, if we use only Fermat reals without incomplete terms. We will 
also see why the name incomplete term seems a better choice than unnecessary 
term. 

4.3 The notion of incomplete term 

Let us start from the usual notations and hypotheses: x = SjLi bj dig^ is the 
decomposition of x, and p G (0, 1]r. In this decomposition, let us consider a 
term dtp r+1 , for r = 1, . . . , N — 1. The power (— ) 1 / p is smooth, because ~ > 1, 
and, with the usual calculations, we can write 

xi = b[ At ph +Ef!) E 1 ' ■ ■ ... ■ b™- 1 dt u(x>l)@p01 , (4.17) 

n=l W 7eN JY ~ 1 7 ' 
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where 



f r /3 2 e/?i n if I ^ N 

k ''~ '~ jmin ( r /3 2 G ft" 1 , i) if | G N 

We have two possibilities to identify the terms, like dt in (dfe + di) 2 , that 
are unnecessary, or, better, incomplete. 

The first one is to say that a term of the type dtp r+1 gives no contribution 
whenever expanding the power x 1 ^, it gives always zero summands, exactly like 
dt in ( d^2 + dt) 2 = dt + ( dt) 2 + 2 • d^2 • dt = dt. Putting it in negative form: 
if, expanding the power x 1 ^, we have that at least one summand, obtained 
from dtp r+1 , is not zero, then the term dtp r+1 gives some contribution, i.e. it 
is necessary. 

A substantial objection against this idea, however, is the following: let us 
suppose that dtp r+1 is necessary, i.e. it gives some contribution. Then, the 
situation described above also includes the possibility that, in the expansion 
of the power x 1 ^, the term dtp r+1 gives, e.g., only one contribution, whereas 
all the other terms involving dtp r+1 give zero. One of our first aims will be to 
prove that, if in the decomposition of x every term gives a contribution, then 
(x^Y = x. The situation can actually be problematic because, following the 
previous extreme example, in x 1 ^ "all the information concerning dtp r+1 " is 
contained in the unique non zero term. For example, 

x = dt 3 + dts_ , dtp T+1 — dt 3 
x 2 = dt 3 + dt 3 + 2 dU ffi 3 = dt 3 + 2 dt. 

2 4 ottJ 2 2 

Can the inverse operation x 2 i— ^ s/x? reconstruct the whole initial information, 
about x, starting only from the unique non zero term 2 dt generated by dtp r+1 ? 



= ckV [l + 2dt lei 



dt 3 ■ I 1 + -2dt 3 - -4dt3 + —8dt 
2 8 5 16 

dt 3 + dt 3 - - dt. 



This counterexample hence gives a negative answer to our question. 

The second possibility of defining a precise notion of incomplete term arises 
from trying to prove the property (.T 1 / p ) P = x starting from a definition based on 



29 



the previous erroneous idea. We will say, intuitively, that dtf} r+1 is incomplete 
whenever expanding the power x 1 ^ gives at least one zero summand. For 
this reason, the term "incomplete" is better than "unnecessary". Putting it 
in negative form: if expanding the power x 1 ^, we have that every summand, 
obtained from dtp T+1 , is not zero, then the term dtp r+1 gives every contribution, 
i.e. it is complete. 

The particular situation of the leading term dtj3 1 is more natural, and is tied 
i i 1 

to the idea that (digjp = dt p p, so that (dtp^" = if and only if j3i < -. 

All this motivates the following 

Definition 34. Under the hypotheses introduced at the beginning of this sec- 

i_ 

tion, we say that dtp r loses information in x? , or that dtp r is incomplete with 
1 

respect to x p , if and only if the following conditions hold: 

1. r = l =► /3i<i 

2. If r > 1, then 

3n=l,... t k x i 3 7 eN 7V - 1 : | 7 |=n, 7r _i ^ , w(x, 7 ) ©pft < 1. 

Let us analyze the condition [2] to see that it corresponds to our intuition: 

• '3n = l,...,k x i 3 7 e N^ -1 : | 7 | = n': looking at (|4.17|) . we can say 
that this part of (^corresponds to "at least one summand in the expansion 
of x 1 ^" 

• ' 7 r-i 7^ 0': "where the term dtp r appears" 

• 'w(x, 7 ) ©p/3i < 1': "is zero". 

Consider the case p = -, where q <G N>o- In this case, the power x 1 ^ 
becomes x q and hence we obtain an equivalent formulation starting from the 
multinomial formula: 

x« = (jrbidtp) = J2 ^f.....bXdt &N _^. 

y=i / lies" 1 " J 

\v\=q 

In fact, we have 

Theorem 35. If p = |, q € N>o, ifit/i 9 < r /?2 Q /^i -1 , i/ien we ftaue i/iai 
dt/3 r+1 loses information in x q 

if and only if 

N 8 

3r) g N w : |r;| = q , r, r+1 ? , ^ < 1. (4.18) 
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Proof: We first compute the term uj(x, 7) © ^ of Definition [Ml in the case 
I7I = q — a, where a € N< g : 



3=1 



N 

1 , Tj-i 



i+E 



A ^ ft eft 

TV 

<? , \ ^ 7j-i 9 , « 



E 



a ^ 2 ft a /v 

where the last equality is due to the hypothesis \j\ = q — a. Therefore 

N 

UJ ( X:1 )( 3 P± = Pl (S mJ^ if| 7 |= ? - . (4.19) 
q a ^ 7,-1 

To prove that (|4.18p is necessary, we start from the hypothesis that there 
exist n = 1, . . . , k„ a = min ( r /32 G ft n , q) = q and there exists 7 £ N N_1 such 
that 

| 7 |=n , 7r7 ^0 , w0,7)©^<1. 

Set 

I g — n if j = 1 



7 l, 



7i _! if./ 2 V 



Then 17 <E N^, = | 7 | + q - n = q, and i] r+1 = 7r ^ 0. Finally, from 1(335]) . 
with a := q — n, we obtain 

which concludes the first part of our proof. 

To prove that (|4. 18|) is sufficient to obtain that d^ r+1 is incomplete, we 
consider 77 as in ((335]) and set 7 := (772, . • . ,Vn) e N^" 1 . Then 

H = \v\ - vi = q - m < 1 = min ( r ft e q) = K± 

1 p 

and 7 r = T] r -\-i 7^ 0. Using (|4.19p with a := r/i < \rj\ = q. we have 
^,7)®^ = ^®© A = ©^<1, 
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which concludes our proof. □ 
The following theorem can be viewed as a validation of our definition of 
incomplete term. 

Theorem 36. Let x 6 D x be a strictly positive infinitesimal, and let p € R>o- 

i_ 

Suppose that in the decomposition of x no term loses information in xp , then 

, /' 



= x. 

Therefore, the power (— ) p is an injection on the set 

| a; G Z?oo | all terms of x are (— )p complete^ U{i/€ *R| °y =/= 0} . 

Proof: If p > 1, setting q := - < 1, the conclusion follows from Theorem l27l 
Therefore, the only interesting case is p < 1. 



If, in the decomposition of x, we have only N = 1 term, then, by hypothesis, 
this term is not incomplete. Taking the negation of Definition [M] for r = 1 , we 
obtain B\ > - and 

/ i\P r IIP ( A \ p 

[xpj = y(bi dt^)" =l6 1 p dt P)3l j =b 1 dtp 1 =x. 

Let us observe that since p(3i > 1, the expression b[ dt p p 1 is a decomposition, 
so that in taking its p-th power, we have applied Definition 1251 We can hence 
suppose N > 1. 

As usual, we refer to (|4.17[) . By hypothesis, every term dtp j is complete, 
which means 

P01 > 1 

and, for every j = 2, . . . , N, n = 1, . . . , k, and 7 £ N^ -1 we must have 

171=71,7^1^0 =*> w(x,~f)®pPx>l. 
This implies that in (|4.17p the greatest infinitesimal is dt p p 1 and every summand 
is not zero. Therefore, to compute ^ifj , we can use Lemma [3"2l 



\xp I 



k / 1 \ 1 



,n l e — ' 7! 

n=l x ' |7|=n 



=: h dt f3l ■ (1 + hf . 
We have to prove that 



h dt Pl ■ (1 + hf = x 

= 61 dtp 1 + b 2 dtp 2 + . . . + bN dtp 
=: bi dtp 1 + k. 

To this end, we use the following 
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Lemma 37. Let h, k 6 Doo and p, b, (3 £ R^o such that p < 1 and (3 > ui(k). 
Define 

k v ^ di 

'■= z_^~r ™a,e/3i 



bdt g ' ^ 6 



where k = 5Z i=1 Oj di Qi zs £/ie decomposition of k. Then 

i 

fc 



l + /i= 1 



od^ 



fodt / 3-(l + ^) p = 6di (3 + fc. 



(4.20) 



Therefore, to obtain our conclusion it suffices to verify the assumption of (I4.20[) . 
that is 



i+ e(!) eS-^-^ 1 ---^- 1 

n=l V 7 kl=n '' 



Indeed, 



b 2 

1 + 7- di/3 2 e^i 



1-yJ 

1 + t- di/3 2 e/3i 

Ol 



"7^- dtfi N Qfi i 



In 



dt 



/l 
P 

, n 



N 



7|=n 
fc /I 



7! \bi 



71 



61 



7W-1 



dt. 



'3=2 7,'-i 



|7| = 



1 + EU) E ^•^"•^••■••<- i d 



as claimed. □ 
Proof of Lemma I37t It suffices to note that the powers (— ) p and (— )~ 

are smooth if applied to invertible Fermat reals. By the elementary transfer 
theorem, they hence have all the usual properties, so that we can write 

n P 



(l + h)> 



1 



bdt 



1 



bdt 







Thus, 



bdt f3 ■ (I + hf = bdt /3 + bdtp 



bdt 13 



bdtp + bdtp ■ y d£ Ql e/3^ 



bdt f 



□ 
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Remark 38. To simplify the notations, let us define 

x is (-)p -complete : -<=>■ Vj = l,...,N: dtp. is complete w.r.t. 

Compl f-J := | a; G \ x is (-)» -complete j U {y G "R | °y 7^ 0} . 
Then the map 

(-) P lco mpl (i) : Compl (£) -H 
is injective. Moreover, the map 

: *R — > *R 

is surjective. However, the power (— ) p doesn't map Compl (^^j onto itself. In 

1 —9 — 

fact, if p = 2 an( l ?/ = ^2, then = y = dt so that every summand in is 

nonzero and hence y G Compl . This y is therefore in the codomain, but it 

is not of the form y = x p for x G Compl (^j , as otherwise y 2 = ~ x y so 

that x = dt and we would also have dt G Compl (2), contradicting ( dt) 2 = 0. 

We can now prove that complete Fermat infinitesimals have very favorable 
properties related to powers. 

Theorem 39. Let x G 13 00 be a strictly positive infinitesimal, and p, q G K>o- 
Then 

x is {—f -complete ==> (x p ) q = x pq . 



Proof: Let us start from the usual formula (|4.11[) applied to x p : 



By hypothesis, x is (— ) * -complete, that is ^ > 1, and for every j = 2, . . . , N, 
n = 1, . . . , fc x . p , and 7 G N^ -1 we have 

H = « , 7J-1 =>■ w(x,7)e^>l. (4.21) 
This property implies that is trivial, in fact, in general 

Va,&GR:a<l =^> (a by 1 = - + \ > - > 1. 

aba 

Therefore 

Vn = 1, . . . , fc^ p V 7 G N^ 1 : | 7 | = n =► wfo 7) > 1, 
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n=l V 7 7 eN N - 1 
7|=n 



■^^l Mffi a. (4-22) 



By (|4.21j) . every summand in (|4.22|) has order greater or equal to 1. We can 
hence apply Lemma l32l obtaining 



(x p ) q = b p 1 q dt i , 
= b pq d* a 

pq 

= b pq dift 



n=l v 7 7G N JV - 1 
|7|=n 



A' 



1 + E^ d< fteft 



3=2 

AT 



1 + E^ d< fte^ 



3=2 



• 6^ 1 dt u ( xn ) 



Here the property [(1 + h) p ] q = (1 + h) pq , for ft 6 -Doo, holds because the base 
is invertible. □ 
The following result does not depend on the notion of complete term, but 
supposes that in the product x ■ y no term becomes zero. 



Theorem 40. Let x, y £ "R>o and p £ R>o, such that 

V* = 1, . . • , N y Vj = 1, . . . ,N X : uj t {y) Wj(x) > 1, 



(4.23) 



where N x and N y are the number of summands in the decompositions of x, y 
respectively. Then 

(x ■ yf = x* ■ y p . 

Let us observe that the hypothesis (|4.23|) reduces to the usual (— ^-comple- 
teness in case x = y, i.e. in case of the property (x 2 ) P — xP ■ x p . 
Proof: We will proceed for x, y £ -Doo, the^proof being analogous if x or y is 
invertible. Let x = $3j=i an d V = Si=i a « dt ai be the decompositions of 

x and y. Then 



x ■ y = ^bjdj dt ai ®p r 

i,3 



(4.24) 
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By hypothesis, every summand in this sum is nonzero. Of course, ot\ ® j3\ is the 
leading term in (|4.24[) , and we can hence apply Lemma [32] to obtain 



{xy) p = a p b p dt ai (Bfi 1 ■ 



1+ E 



\ 



Mi 



at 



a,efte(aieft) 



On the other hand 



N 



3=2 



f = a p dt sx - l + V — dt ajeai ] . 

p V / 



The equality (|4.25|) can also be written as 



(4.25) 



(4.26) 
(4.27) 



(xy) p = a P b P dt ai ®p 1 - ( 1 + ^ dt aiefte(ctlffii8l ) 



Mi 
6iai 

Mi 
Mi 



EM* ,, 
7— ar Qi©^ie(ai©/3i) 



■ E 7TT" d *"i®fte(ai©/3i) 



i>2 

which is the same result obtained from multiplying (|4.26[) and (|4.27p . □ 

4.4 Roots are not smooth 

Here we prove that the power function 

(—) p : *R>o — > *M , P<1 

is not (non standard) smooth on D x \ {0} (for the notion of non standard 
smoothness, see |1Q[ [5]). This is a naturally expected result, because the corre- 
sponding derivative should be 

dx w /I 1 "? 

and, intuitively, . fL p is an infinite, whereas in *R we obviously do not have 
infinities. Therefore, the theory of Fermat reals should be sufficiently complete 
to prove that the power function (— ) p is not smooth at h E Doo \ {0}. Indeed, 
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from the generalized Taylor formula (see Theorem 12.1.3 and Definition 12.2.7 
in [9]), if we assume that (— ) p is smooth, we would have 

+ oo ,j 

3m £ *K N : Vfc £ D M : (h + kf = ^ — • mj. (4.28) 

3=0 :h 

For the sake of completeness, we recall that this sequence m = {rrij)- e¥i is 
unique up to first order infinitesimals, i.e. if m = (mj)j £K verifies (|4.28p . then 
nij =i fhj for every j £ N. 

Set s := -sgn(/i) e { + 1,-1}, and take k = sdt in (|4.28p . We have that 
h + sdt 7^ because sgn(h) ^ sgn(sdi), and taking into account that mo = 
P = 0, we obtain 

(h + sdtf =m + sdt-mi = sm 1 dt G h = {x E "M| u{x) < 1}. 

Now, the order of (h + s dt) p is max (^( /i )- 1 ) = f£^2 because h + sdt =^ and 

since always uj(h) > 1. Therefore, we would have < 1, i.e. p > uj{h) > 1, 
whereas we have supposed p < 1. 

4.5 An application to the infinitesimal Taylor formula with 
fractional derivatives 

Using powers hP of infinitesimals h £ -Doo, we can prove an infinitesimal Taylor 
formula with fractional derivatives in a straightforward manner. This further 
underlines the ease of translating classical results using the infinitesimal lan- 
guage of the ring of Fermat reals. Frequently, these translations are really 
faithful to the informal use sometimes appearing in applications. Let us note 
that the same translations are not so easily performed in algebraic models of 
infinitesimals, like in Synthetic Differential Geometry (see, e.g., [21], [15]) or in 
Levi-Civita fields ([26] [25]) or Weil functors ([Ml [18]). 

We start with some definitions and a theorem, taken from [23 . 

Definition 41. If a £ ffi, we will denote with C Q (IR>o,IR) the set of all the 
functions / : R>o — > K that can be written as 

f(x) = x p ■ h(x) VxGM> , 

for some p > a and some continuous function /i £ C°(R>o,R). Moreover, for 
every m £ N>o U {oo} we also set 

C'(R>o,K) := {/ £ «>o — > R| / M € C Q (M >0 ,M)} . 

Secondly, we define the Riemann-Liouville integral operator of order a > 
with a > 0. 

Definition 42. Let a £ R >0 , a £ M> and / £ C Q (M >0 ,M), then 
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J°J(x) := f{x) 

Jaf(x) := / ix- t) a -\f{t) dt Vx g K >Q . (4.29) 



r(«) 7a 

Here T denotes the gamma function. To derive the fractional Taylor formula, 
we need the Caputo fractional derivative. 

Definition 43. Let a g R >0 , a g R> , and / g C m 1 (M >0 ,K). For simplicity of 
notations, let m := r a~ l be the ceiling of a. Then 

£>«/ : R> Q — ► R (4.30) 

D%f(x):=J™- a fW(x) Vx>a. (4.31) 

Finally, we set 

From (|4.30D and (|4.29|) we therefore have 

= p7 1 x A* - tr-^f^Ht) dt, 

T(m-a) J a 

where m = r a~ l . 

The non-infinitesimal version of the generalized Taylor formula with frac- 
tional derivatives is the following. For its proof, see |23| . 

Theorem 44. Let a, a, b g R, and n g N, with < a < b and < a < 1. 

Consider a continuous function f g Co([a, &],M) smc/i i/ia£ 

D*' a / gC ([a,6],R) Vfc = 0,...,n + 1. 
T/ien /or every x g (a, 6] £/iere exists £ g [a, sc] suc/i £/ia£ 



- L f(iaTT) a /(a) + F((n + l)a + l) ' ( * " ° } 



In our framework we are able to prove a corresponding infinitesimal Taylor 
formula for the following class of smooth functions: 

Definition 45. In the hypothesis of the previous theorem, we set 

C~([a,&],R) := {/ g C°°Qa,b],R) \ D k a > a f g C (M],R) Vfc g N} . 

Finally, we can state the main result of this section: 

Theorem 46. Let a, a! , a, b € R, and n g N, with < a' < a < b and 
< a < 1. Consider a smooth function f € C£°([a, &],R), i/ien 

™ Lie* 

V/i g : /(a + h) = ^ -— __Dj«/(a). 
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Proof: Let h = Xw=i dtp } be the decomposition of the infinitesimal h £ 
D( n +i)a- By Definition 1251 we have 



h ia = b\ a dt^ ■ ( 1 + J2 £ difteft J Vi = 0, . . . ,n. 

N r 

This means, using an innocuous abuse of language, that ht := Ylj=i bjffi and 
K\ a := b\ a t%T . [ 1 + ^ ^t 1 ^ ] Vi > 

are little-oh polynomials representing the Fermat real h and h ta respectively, 
and 

(h t y a = h] a vi > o. 

For t > sufficiently small, we have a' < a + h t < b, and we can apply Theorem 
H4"lat the point a, obtaining 



2 — 

(4.32) 

Now, /i G r>(„+i)a.-i so that < (71 + l)a and /7,(™ +1 ) Q = 0, that is 

,(n+l)a 

U m 13 = 0. (4.33) 

Considering that _D™ +1 ' Q / is continuous on [a, b] and that £ t £ [a, 6], from 
([435)1 and P^3"2"j) we obtain 

n uicx 

f{a + h t )=Y jVF ± -T,DTf{x)+o{t) asi^0+, 
^r(za + l) 

yielding the claim. □ 



5 Computer implementation 

The definition of the ring of Fermat reals is highly constructive. Therefore, 
using object oriented programming, it is not hard to write a computer code 
corresponding to *R. We realized a first version of this software using Matlab 
R2010b. 

The constructor of a Fermat real is x=FermatReal(s,w,r), where s is the 
n + l double vector of standard parts (s(l) is the standard part °a;) and w is the 
double vector of orders (w(l) is the order uj(x) if x € *R\R, otherwise w=[] is 
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the empty vector). The last input r is a logical variable and assumes the value 
true if we want that the display of the number x is realized using the Matlab 
rats function for both its standard parts and orders. In this way, the number 
will be displayed using continued fraction approximations and therefore, in many 
cases, the calculations will be exact. These inputs are the basic methods of every 
Fermat real, and can be accessed using the subsref , and subsasgn, notations 
x.stdParts, x. orders, x.rats. The function w=orders(x) gives exactly the 
double vector x . orders if x <G *R \ R and otherwise. 

The function dt(a), where a is a double, constructs the Fermat real dt a . 
Because we have overloaded all the algebraic operations, like x+y, x*y, x-y, -x, 
x==y, x~=y, x<y, x<=y, x~y, we can define a Fermat real e.g. using an expression 
of the form x=2+3*dt(2)-l/3*dt(l), which corresponds to x=FermatReal ( [2 
3 -1/3] , [2 1] ,true). 

We have also realized the function y=decomposition(x) , which gives the 
decomposition of the Fermat real x, abs(x), log(x), exp(x), isreal(x), is- 
inf initesimal (x) , isinvertible (x) . 

The logical function v=eqUpTo(k,x,y) corresponds to x y. 

The ratio x/y (see Theorem I22[) has been implemented for x and y infinites- 
imals and y~=0, or in case y is invertible. Finally, the function y=ext(f ,x), 
corresponds to ' f(x) and has been realized using the evaluation of the symbolic 
Taylor formula of the inline function f . 

The functions dF and dOmega correspond, respectively, to the Fermat and 
the omega distance, while x~p, sqrt(x) and nthroot (x,n) have been realized 
both for x infinitesimal or invertible using the formulas we have derived in the 
present work. 

Using these tools, we can easily find, e.g., that 



sin(Vdt 3 + 2dt 2 ) , , 2 , 1096 , 1234 , 
cos(^CT) " * + * - 3* + m? d 't + HT* 

This corresponds to the following Matlab code: 
>> x=sqrt(dt(3)+2*dt(2)) 
x = 

dt_6 + dt_3 - l/2*dt_2 + l/2*dt_3/2 - 5/8*dt_6/5 
>> y=nthroot(-4*dt(l) ,3) 

y = 

-1008/635*dt_3 
>> g=inline ( ' cos (y) ' ) 
g = 

Inline function: g(y) = cos(y) 
>> f =inline ( ' sin(x) ' ) 
f = 

Inline function: f(x) = sin(x) 
>> decomposition(ext (f ,x)/ext (g,y) ) 
ans = 

dt_6 + dt_3 - 2/3*dt_2 + 1096/2787*dt_6/5 + 1234/913*dt 
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Up to now, this code has been written only to show concretely the possibil- 
ities of the ring *R. On the other hand, it is clear that it is possible to write 
it with a more specific aim. For example, as in case of the Levi-Civita field 
([31[2H|) possible applications of a specifically rewritten code include automatic 
differentiation theory. Let us note that, even if the theory of Fermat reals ap- 
plies to smooth functions, a full treatment of right and left sided derivatives is 
possible so that the theory can be applied consistently also to piecewise 

smooth functions. Finally, the use of nilpotent elements permits to fully justify 
that every derivative estimation of a computer function f[25|) reduces to a finite 
number of algebraic calculations. 

The Matlab source code is freely available under open-source licence, and 
can be requested from the authors of the present article. 

6 Conclusions 

Usually, it is common to study extended structures, like the ring of Fermat 
reals *R, using suitable extensions of well established notions. For example, it 
is more natural to search for metrics of the form d : *R x *R — > *R than for 
standard metrics on the set *R. We have shown that it is possible, and also 
very natural, to define standard topological structures on the ring *R having 
very favorable relationships with various aspects of differential calculus on *R. 
This allows a better dialog with mathematicians not already familiar with the 
theory of Fermat reals and underlines that this ring is not "non standard". 

Moreover, with the present work, we are continuing our program to define a 
meaningful and powerful ring with infinitesimals using only very well behaved 
representative functions for new numbers. If one thinks at non standard analysis 
or Colombeau's ring of generalized numbers, it becomes clear that this quest is 
nontrivial. As a consequence, we have been able to characterize ideals of the 
ring "R in a very simple and descriptive way. 

Finally, we have proved that nilpotent elements and arbitrary roots can 
coexist very well, even if this seems impossible at a first glance. This is a very 
important step toward the idea of using nilpotent infinitesimals for stochastic 
calculus. For example, based on a very helpful discussion with N. Blagowest 
(Department of Physics, K. Preslawki University, Bulgaria) we may call Ho 
process any (deterministic) function x : *R — > *R such that 

yt e R V7i e Di : x(t + h) = x(t) + v [t, x(t)] ■ h + \Vh. 

In this approach, the deep mathematical problem is that there doesn't exist a 
non trivial smooth function that verifies such a definition. Of course, we need 
continuous but nowhere differentiable functions and hence, we need to extend 
the ring *R by suitable infinities. Indeed, this is indispensable if we want to 
consider the derivatives of the function x. Therefore, the new problem to face 
becomes: can infinities coexist as reciprocals of nilpotent infinitesimals? This 
question will be the subject of future work. 
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